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Abstract 

A holographic model of QCD at finite temperature is proposed and ther¬ 
mal properties for light mesons are examined in terms of this model. We find 
interesting temperature dependences of masses and decay constants for several 
mesons, and Nambu-Goldstone theorem for the chiral symmetry breaking and 
generalized Gell-Mann-Oakes-Renner relation are assured at finite temperature. 
Furthermore, we can see in the chiral limit that the speed of the pion in the 
thermal medium decreases to zero when the temperature approaches to the crit¬ 
ical point of quark confinement and deconfinement transition, where the chiral 
symmetry is restored. 
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1. Recently, from the gravity/gauge correspondence, probe brane approach to the 
gauge theory with flavor quarks has been largely developed in terms of the system of 
D p /D p+4 branes [Tj, 2]- bi these, flavor quarks have been introduced by embedding 
D p .|_ 4 branes as probe in the lOd background, which is deformed to the anti-de Sitter 
space (AdS) by the stacked D p branes. Then several important dynamical properties 
of QCD as chiral symmetry breaking and quark confinement have been explained, and 
furthermore, the flavored-meson masses have been calculated. 

Inspired by these works, phenomenological, five-dimensional holographic models 
have been proposed to explain more quantitatively the physical properties of light 
mesons in a simpler setting [2 El U HUE!- Up to now, all these works of this framework 
are restricted to the case of zero temperature (T = 0). Meanwhile, it is interesting to 
extend this simple framework to the finite temperature scheme in order to compare the 
thermal properties of the quantities estimated in this framework and the recent exper¬ 
iments jSJ. Up to now, such thermal properties of light mesons have been estimated 
by the lattice, chiral perturbations and sigma models. So, we extend here the model 
proposed at T = 0 in 32 m to the finite-temperature version in order to give the first 
holographic approach in this category. 

bi 32 S3 , the gravity dual of the gauge theory with confinement at T = 0 is ex¬ 
pressed by the 5d AdS with an infrared (IR) cutoff at an appropriate position of the 
fifth coordinate. This IR cutoff is essential to form a gravity dual of a gauge theory 
in the confinement phase, and the position of IR cutoff is very sensitive to the me¬ 
son mass. In this sense, this point is fixed rigidly. The finite-temperature version is 
analogously obtained by replacing the AdS background with the AdS- Schwarzschild, 
and the horizon, which is related to the temperature of the gauge theory, in the AdS- 
Schwarzschild is set inside the IR cutoff at low temperature, then the horizon is hidden 
by the IR cutoff in the 5d bulk. 

It is well-known that the 5d AdS-Schwarzschild background corresponds to the 
finite temperature version of CFT with deconfinement phase when the IR cutoff is 
removed. This situation is realized in the present case when the IR cutoff is put inside 
of the horizon. In this case, the IR cutoff is covered by the horizon, then the effect 
of the IR cutoff disappears. This setting corresponds to the high temperature phase 
of our model. On the other hand, at low temperature phase, the cutoff is set outside, 
the fields in the bulk end at this IR cutoff and appropriate boundary conditions for 
the fields are imposed at this end point. As a result, we find infinite series of discrete 
meson spectra. This situation is completely parallel to the AdS case, but the system 
is set at finite temperature in the present case. 

The temperature is understood from the finiteness of the Euclidean time in the 
AdS-Schwarzschild configuration in order to evade the conical singularity which may 
arise at the horizon. This constraint for the Euclidean time should be preserved, even 
if the IR cutoff is introduced, in order to connect the manifold smoothly up to the 
critical temperature of quark confinement and deconknement (T c ). So we can perform 
the analysis for the properties of the light mesons in the region of 0 < T < T C1 by this 
configuration. As for the quark part, we use the effective 5d action (SAeson) introduced 
in 1203 for the sake of the simplicity. 

2. The bulk background, which is dual to the finite temperature 4d QCD, is formed 
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by introducing the IR cutoff HUE) z m into the 5d AdS-Schwarzschild solution, 

d 4 = i (^-f 2 {z)dt 2 + (dx l ) 2 + , f( z ) = 1 - (y^) 4 • (!) 

where the radius of AdSs is taken as unit. As is well-known, the imaginary time is 
restricted as 0 < it = t < ttzt, to avoid the conical singularity. Then the temperature 
T is given as T = 1/(7 tzt)- As for the infrared cutoff z m , we consider the case of 
z m < zt in order to see the confinement phase at finite temperature. I 11 this case, fields 
are restricted to the region 0 < z < z m . 

The condition 0 < it = t < 7 tzt is not altered even if the cutoff is introduced. Then 
the temperature is retained as given above. When the cutoff has been introduced, 
however, the possible conical point is removed for z m < zt, so it may be possible to 
change the temperature in order to find a smooth manifold for z < z m . But here we 
fix the temperature as T = 1/(7 tzt) in order to extend smoothly the region of z t up 
to Zt = z m . The confinement is thus assured since we have infinite series of discrete 
meson spectrum, quark and anti-quark bound states, due to the boundary at z m . 

On the other hand, for the case of z m > Zt, the situation is the same with the 
finite temperature CFT. In this case, the cutoff has no meaning since z is restricted to 
Zt > z > 0. Then we find the deconfinement phase with a few number of possible meson 
states, whose masses should be below the double quark-mass. This point is assured in 
|2]. This implies that the critical temperature T c of confinement-deconfinement is set 
in the present case at T c = 1/(7 Tz m ) which is estimated here to be ~ 102 MeV from 
meson spectra. This value is not so wrong since our model is not so strict. 

3. For the sake of the simplicity, according to [HI IE], we start with the following 
5d meson-action under the background ©. This action is considered as the quark 
part represented by the probe brane in the string theory. And the fields on the brane 
represent meson states regarded as the bound states of quark and anti-quark. The 
fields considered here are the gauge fields, Lm and Rm, and a scalar held <f> whose 
vacuum expectation value (VEV) is connected to chiral symmetry breaking. Here, we 
consider Nf = 2 flavors, and $ transforms as a ( 2 l, 2r). The action is 
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where the covariant derivative is defined as D M Q = 8m& + iLm^ ~ i&Rm , 9 is the 
determinant of the metric, and Lm = L a M r a for the Pauli matrices r a and similarly 
for other fields. g$ is the 5d gauge coupling. We define <f> = S e mara and \v(z) = (S), 
where S corresponds to a real scalar and 7r to a real pseudoscalar (S S and 7 r —> — 7r 
under L R). They transform as 1 + 3 under SU(2)y. 

3-1. The 5d mass of the scalar is set as M| = —3 to consider the bulk held which 
corresponds in the gauge theory side to an operator with the conformal dimension 
A = 3. Then the equation of motion for v is given as, 
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And we obtain 
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Here two integral constants, m q and c, are identified with the quark mass (explicit 
breaking of the chiral symmetry) and the chiral condensate (spontaneous breaking of 
chiral symmetry in the chiral limit), respectively. Actually, we get the approximate 
form, v(z) ~ m q z + cz 3 near z = 0, which is consistent with the above statement. 

We notice here that v(z) varies with temperature T or zt, which is specific to our 
model, but the parameters m q and c are independent of T. In other words, the gauge 
theory considered here is characterized by these three parameters. However, the chiral 
condensate c should be determined dynamically and it should depend on T. This 
point is the defect of the present model, and we must be careful about the results 
of our calculations for the quantities, which depend on this ’’parameter” c, when we 
discuss the T dependence of them. 


The fluctuation of S, which is defined as S = v(z)/2 + a, can be observed as a 
singlet meson state (a). Here and hereafter we consider the static mode, d t (f) = 0, for 
any field / in order to derive the mass in a simple way. So the invariant mass, or pole 
mass, is defined here as —<9 2 / = m 2 (j). Then the equation for a is given by adding this 
4d mass term to Eq. @ as 
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It should be noticed that this equation is independent of v(z). In this sense, the 
model given here is insufficient since all the fluctuations on the probe brane would be 
influenced by the shape of the embedded brane. But the improvement of this defect is 
postponed to the future. 

The discrete mass-spectrum is obtained by solving this equation with the boundary 
conditions, cr(z)\ Z0 = d z cr(z)\ Zm = 0, where Zo is the UV cutoff which is taken to zero 
after all. The mass depends on z m and T. By using the value of z m determined by the 
p meson mass at T = 0 as 1 /z m = .323 GeV |6J, the T-dependence of the a meson mass 
is estimated. The numerical results are shown in the FiglU We can observe that the 
mass decreases with increasing T and approaches to zero for T —*■ T c . This behavior 
is easily understood from Eq. ©• The mass term, m 2 // 4 , in the equation is largely 
enhanced near z = zt, where / —> 0, so m should be very small when z approaches 
to Zt■ This situation is realized for zt —*► z m or T —> T c . This point is therefore the 
most prominent check point of this model. The experiments to assure this property 
are acquired. 


3-2. The gauge bosons are separated to the vector and the axial vector bosons Vm 
and Am, and are defined as L M = Vm + A m and Rm = Vm — A M , respectively. 

First of all, we consider the vector mesons. At finite temperature, it is impossible 
to write the equation of motion covariantly, so we consider the linearlized equation for 
the spatial component V) here. And it is (employing V z = 0 gauge) given as 
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This equation gives us the discrete eigenvalues m 2 = m 2 . In order to see the decay 
constants, it is convenient to expand Vi{x,z) as Vi(x,z) = \ x )^n( z ) f° r eac h 

mass eigen-state. The spacial components of the wave-function for each mode are 
normalized as 


dz 


IZO 


zp 


= 1, 


(7) 


where we notice the factor 1/f 2 (z) in this integrand. This is an important factor in 
determining the T dependence of the decay constant near T = T c . Actually, we find 
rapid decreasing of F p near T = T c due to this factor as shown in FigH] 

The mass of vector meson is given by the boundary condition Vi(z 0 ) = d z Vi(z m ) = 0 
similarly to the case of the sigma meson. Note that Eq.© is also independent of v(z), 
then the masses of vector meson depend only on T, and the T-dependence of the lowest 
one, which is identified with p meson, is shown in FigjT] 


m [MeV] 



VF [MeV] 


T [MeV] 



T [MeV] 


Fig. 1: (A) The curves show the T dependence of m v , m p , m a and m ai from the 

bottom. In (B), the solid curves denote the T dependence of F^. and yjF ai 

from the bottom. The dotted curve denoting F% is constant with high accuracy at 


T < 100 MeV. 


We can see that the mass of p meson vanishes or becomes very small near T = T c . 
This is consistent with the expectation obtained from the equation © as in the case 
of the a-meson. 


3-3. The linearized equations of motion for the axial vector meson A p and the 
pion 7r are obtained by decomposing A p into the transverse and the longitudinal part, 
A p = A p ± + d p (p] for simplicity, flavor index is neglected. Using A z = 0 gauge, one 
then one obtain 
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where the coupling constant g§ is determined from the vector current two-point function 
at T = 0 0. Here m a and are the pole masses dehned as —d 2 ip = m 2 (p and 
—<9 2 7r = m 2 7r. These equations are solved numerically under the boundary conditions, 
A^-l( z o) = = 0 and (p(z Q ) = d z p(z m ) = i r(zo) = 0. 


The decay constants of the axial mesons and the pion are calculated from the wave 
functions as 0 


F 2 = — 
9l 


d 2 h£ 


dz 2 

ZO- 


1 2 


(n % 0), 


i 9,41 

(F * > = ~n — 


ZO 


( 12 ) 


where Ai±(x, z) = ’}^ n ai L \x)h^ L {z) and h^(z) is normalized as h^(z) given in (0. 
Furthermore, A^f (A-^) is the zero-mode solution to Eq. 0 (Eq. 0 with rn 2 = 0), 
satisfying Aq 0 - ± (z 0 ) = 1 and d z Ajj > j ± (z m ) = 0, and F* ,s are the timelike and spatial 
components of the pion decay constant. Obviously, F \and F£ are different from each 
other at finite T, while F* — F* at T — 0. 

The masses and the decay constants of axial-vector mesons and pion depend on 
four parameters m q , c, z m and T through v(z) and f(z ) in Eqs. (l8l)- (fHl) . while those 
of vector mesons do only on z m and T. For the consistency between the vector and 
axial-vector meson sectors, here we take the same z m as that determined in the vector 
meson sector; namely l/z m — 0.323 GeV. Parameters, m q and c, are determined to 
reproduce the experimental values, rh n and F n , of m n and F n at T = 0; the resultant 
values are m q = 2.26 MeV and c = (0.333 GeV) 3 . This parameter set reproduces 
masses and decay constants of 7r, cr, p, cq mesons at T = 0 within ~ 10% error M. 
Furthermore, when T = 0, the present model satisfies the Gell-Mann-Oakes-Renner 
(GOR) relation fn\F 2 = 2 m q c 0 6]. 

In principle, the chiral condensate c depends on T, but it is known by lattice QCD 0 
and the chiral perturbation theory to three loops ra that c little depends on T in the 
region of T < 100 MeV. So we focus our discussion on the temperature region for a 
while and then neglect the weak T dependence of c. 

Figure |U shows that the T dependence of and Ff is almost constant in the region 
T < 100 MeV; more precisely, is slightly reduced and F is slightly enhanced there. 
This is consistent with the result of the chiral perturbation theory to two loops HU. 
Furthermore, m a is largely reduced as T increases, as shown in Fig. |U This is also 
consistent with the result of lattice OCD IQ; note that the pole mass is evaluated in 
the lattice calculation. Thus, the present model qualitatively simulates T dependence 
of real QCD in the region T < 100 MeV. 

In Fig. HJ the T dependence of m a is relatively weaker than that of m p . The 
difference stems from a term —(g 5 v) 2 /(zf) 2 which is present in Eq. 0 for the axial 
vector meson but not in the corresponding equation 0 for the vector meson. The 
term suppresses partially the enhancement of the mass term m 2 // 4 near T C) since the 
two terms have opposite signs to each other. In Fig. we can also see a tendency 
for m n to become small as T reaches T c . When T « T c , the second term of Eq. m 
is enhanced by the factor / 2 compared with the first term. This might indicate that 
tp(z) ~ Inserting this relation into Eq. (fTTl) . we can speculate that m n cx vf 

near T c . 
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m 2 [MeV 2 ] 



[MeV] 


Fig. 2: The m q dependence of m); at T = 100 MeV. The solid line corresponds to 
m%(T) and the closed circles to 2 m q c/F* 2 (T). 


3-4. In real QCD, a generalized GOR relation, m%(T)F% 2 (T) = 2 m q c, is satisfied 
for finite T |[THl lllj . Within the framework of the holographic QCD, a necessary 
condition for the GOR relation to be satisfied is that Eq. © agrees with Eq. hud 
with n(z) = 0 PJ. This condition is satisfied for any T in the present model. Thus, the 
present model is expected to reproduce the generalized GOR relation properly. 

Whether the generalized GOR relation is realized in the present model is tested 
through m n (T) and F*(T) calculated numerially. Figure. El shows the m q dependence 
of m^(T) and 2 m q c/F* 2 (T) by the solid line and the closed circles, respectively. Here 
the case of T = 100 MeV is taken as an example. The m^(T) shown by the solid line 
tends to zero as m q deceases, as a reflection of the Nambu-Goldstone theorem at finite 
T. Comparing the solid line and the closed circles, one can see that the generalized 
GOR relation is well satisfied in the present model. 


3-5. The pion velocity v n in the thermal midium can be estimated from F*{T) and 
F‘(T) as a, 2 = F;(T)/F‘(T) d3|H,. The pion velocity v n thus evaluated decreases 
quite slowly from 1 to 0.97 as T increases from 0 to 100 MeV; here we assume c does 
not depend on T. Now the value of v n is estimated in the limit of T —> T c by solving 
the Eqs. © and (HDD at the chiral limit rri q = 0. At T c , we expect the restoration of 
chiral symmetry, so we set as v(z) = 0 in the equations and the solutions are given as 


Ai — a\ + a2 


tanh ( z/ztY 

24 


Aq — b\ + 


(13) 


The constants 01,2 and W 2 are determined by the boundary conditions, A ij0 |z=o = 1 
and d z A i>0 \ z = Zm = e, and we take the limit e —> 0 after obtaining v n . Thus we get 


vl = ^1 - (T)“ (X i 1/2 , (14) 

where t = (T c — T)/T c . This indicates v n (T c ) = 0. This is consistent with the result 
obtained by Son-Stephanov M- The corresponding critical exponent is v — 1/2 in our 
case. 
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4. The model we present is the first holographic model to describe QCD in the 
region of T <T C nonperturbatively. In this model, we take the 5d AdS-Schwarzschild 
solution as the gravity background, and cut off the extra dimension at an appropriate 
infrared point z m in order to introduce the confinement in a gauge theory dual to the 
gravity. The present model then describes the confinement phase in the region T < T c . 

The present model makes it possible to evaluate T dependence of meson masses and 
decay constants nonperturbatively. In the region T < 100 MeV, the pion mass and 
the pion decay constant little depend on T. This is consistent with the result of the 
chiral perturbation theory to two loops 0 In contrast, the cr meson mass is largely 
reduced as T increases. This is consistent with the result of lattice calculation P2|, 
Furthermore, the Nambu-Goldstone theorem and the generalized GOR relation are 
assured at any finite T. Thus, the present model reproduces the qualitative properties 
of QCD in the region T < T c . We then conclude that the present model is reliable 
enough to make qualitative predictions on QCD in the region T < T c . 

Our prediction of the pion velocity is v w (T c ) = 0. This agrees with the result of Son- 
Stephanov m based on scaling and universality. Thus, the present result also supports 
the statement that the measured pion velocity 0.65 ra. deduced from the pion spectra 
observed by STAR 0 at RHIC, would be a signal of QCD phase transition. 

The present model has three parameters; quark mass m q , the chiral condensate 
c and the cutoff z m . In principle, c and z m should be determined dynamically and 
be functions of T. The T dependence of these quantities can not be obtained by the 
present model. The lattice calculation shows that the pion mass somewhat increases 
near T c . If it is true, this would imply that z m becomes smaller there, since smaller 
z m yields larger pion mass. Furthermore, it is well known that c has a strong T 
dependence near T c . This is essential for properties of light mesons near T c . Thus, it 
is quite important to improve the present model so that c and z m can be determined 
dynamically. Furthermore, the vector meson mass and decay constant do not depend 
on c in the present model. This is also an important issue to be solved in future. 

K. G thanks Y. Nakano for enlightenment on finite temperature field theory. This 
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of the Ministry of Education, Science, Sports, and Culture of Japan. 
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